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The reliability and service life of power cables is closely related to the cable 
ampacity and temperature rise in the power cable. In a conductor carries AC 
current, complex processes may appear, which directly affect the temperature 
of the conductor surface. So, to keep a conductor in a good state, it is necessary 
to maintain the conductor temperature in a acceptable value. In this paper, a 
procedure for measuring the temperature of conductor surface and the 
corresponding numerical processing of measurement results has been 
presented. The measurement of the temperature probe characteristics and the 
temperature measurement on the surface of the conductor, both required the 
use of certain numerical methods, such as interpolation and fitting of the 
measured values in time diagrams. The procedure was applied to three copper 
conductors with different cross section area and one aluminum conductor and 
the final results are presented graphically, in the form of time diagrams. 
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1. INTRODUCTION 

In addition to the development of intelligent networks, sustainable operation of high-voltage cables, 
as an important career in electricity transmission, contributes to guaranteeing the quality of energy resources 
and improving the safety and security of transmission lines [1]. The temperature of the alternating voltage 
conductors is one of the important parameter that directly determine their endurance [2]. Without a detailed 
study of the current working environment for high voltage cables, raising the conductor temperature above the 
allowable value will accelerate the aging of the cable insulation, increase the leakage current and ultimately 
lead to insulation failure [3]. Therefore, the temperature monitoring of cable conductors can estimate the 
carrying capacity of cables in real time and guarantee their safe and efficient operation [4]. However, in the 
case of time-varying currents, there is a surface effect in the conductor and an uneven distribution of current 
across the cross section of the conductor, so the problem is further complicated [5]. The specific resistance of 
the conductor changes with the change of temperature, which also affects the current distribution. The whole 
phenomenon is not linear, so a very complex mathematical apparatus is needed to determine the temperature 
distribution on the surface of the conductor, caused by the complex shape of the Joule effect [6]. 

During recent years, the development of online monitoring technology for the surface temperature of 
high voltage cables has become a standard for knowing the conductor temperature in high voltage cables [7]. 
There are two methods for monitoring the temperature of the transmission cable, which are contactless and 
contact temperature measurements. The infrared measurement technology method takes faster response time 


Journal homepage: http://beei.org 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 O 2511 


and requires simple equipment structure, but at the same time it is affected by infrared electromagnetic waves 
and environmental factors, which leads to poor measurement accuracy [8]. 

In the case of contact temperature sensors, thermocouple sensors are used to perform multipoint 
measurements. Compared to infrared thermometry, the measurement is flexible and allows temperature to be 
measured in several locations, but it is more difficult to maintain a safe electrical distance between sensors and 
equipment with higher voltage levels. Therefore, neither of the above two methods is considered the most 
suitable option for long-distance continuous monitoring of the temperature of cable conductors in complex 
laying environments [9]. These methods are used to find out the surface temperature of the conductors of high- 
voltage cables, but direct measurement of the conductor temperature is difficult to be achieved. 

Currently, the main methods for calculating conductor temperature used in power cables include 
equivalent thermal circuit models [10], artificial intelligence algorithms [11], in addition to finite element 
analysis [12]. Also the standard (IEC 60287) is the standard used for the equivalent thermal circuit method, 
which gives a complete analysis method and an integrated formula for calculating the load capacity [13], but 
it does not take into account the effect of instantaneous changes in the external environment and other factors 
affecting the correction factor and therefore it cannot solve physical problems in an effective way, such as 
convection and radiation in addition to the coupling of heat transfer [14]. This leads to calculation errors that 
cause the final results to deviate from the actual values in complex and multi-loop environments therefore, it 
is not suitable for monitoring the surface temperature of conductors in real time [15]. 

Also, artificial intelligence algorithms have been applied for the purpose of predicting the temperature 
range of conductors by using support vector machines (SVM) [16], the trans-conductor temperature model was 
developed and the particle swarm optimization (PSO) algorithm was introduced in order to improve the 
parameters of the network model. The surface temperature of the conductor used in the measured cable and the 
load current were also used as basic model inputs to obtain the dynamic temperature of the conductor [17], but 
if this method is used, a large amount of sample data is required for the proposed model to calculate the 
conductor temperature [18]. In this paper, the numerical processing of the temperature measured values on 
conductor surface have been conducted by using MATLAB software package on a three samples of copper 
conductors with different cross-sectional areas, and on one sample of aluminum conductor. 


2. METHOD 

In the case of an AC current flowing through a conductor, uneven distribution of current across the 
cross section of the conductor will occur [19]. When the conductor has a cylindrical shape, it is most convenient 
to solve the problem in a cylindrical coordinate system, so that the axis of the conductor is placed in the 
direction z of the coordinate of the adopted system [20]. So, the vector of electric field and the current density 
vector have only z component, and the magnetic induction vector has only ọ component [21]. These quantities 
depend only on time (t) and distances from the conductor axis (r) and do not depend on the coordinates @ and 
z as shown in [22]: 


E, = 177. Ez (1) 
J=tJz (2) 
H = Tg. Ho (3) 


where: E, is the vector of electric field strength in the conductor, J is the vector of current density and H is the 
vector of magnetic field strength. In the case when the current changes over time according to the simple 
periodic law, a complex notation can be used to solve the problem. In this case, Maxwell's first two equations 
in complex scalar form are used as shown in [23]: 


-F= -poH (4) 
(rH) =J (5) 


where y is the permeability of the medium and œ is the angular frequency of the field source. When the value 
of H in (3) is substituting in (4), a partial differential equation for the spatial distribution of the complex 
intensity of the current density vector can be obtained as shown in [24]: 


dJ 
+p = Jpwo] (6) 


rdr 
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where o is the specific conductivity of the conducting material. The Bessel differential equation of zero order 
can be obtained as shown in: 


d? d 
ae ae (7) 


du? udu + 


U= J's (8) 


A (9) 
K = ono (10) 


now, the distribution of the intensity of the current density vector over the conductor cross section J(0) can be 


expressed through the Bessel function as shown in: 
J(s) = J(0)ber(s) + jbei(s)] (11) 


where ber(s) and ber(s) are the real and imaginary part of the Bessel function which represent the specific 
resistance of the conductor changes, more or less, with a certain temperature (9) for smaller and larger 
temperature changes as shown in [24]: 


Ps = Po(1 + ad) (12) 
Po = Pol + a + Bo +yd*) (13) 


where a, B and y are the temperature coefficients of the conductor material, pois the specific resistance of the 
material at 0 C. Expressions for changes in specific resistance due to changes in temperature in AC current can 
be written as shown in [24]: 


Ps = Po(1 + a208 — 20))(1 + Ys + Yp) (14) 
Ps = KıK2P20[1 + a29(8 — 20)] (15) 


where k;, ys are influence coefficients of the surface effect on the specific resistance of the material, and k2 is 
the losses influence in ferromagnetic material. So, the volume density of Joule losses at each point of the 
conductor (dv) will be: 


dPy 
dv 


Joule losses P; can be determined by integrating expression (16) over the total volume of the observed 
conductor as shown in: 


P, = Jf -av = Jff J : E.dv = {ff p: J? dv = fff o: E? -dv (17) 


the energy converted from to heat during the time interval t, will be: 


dP 


W, = fy Sf E -dv - dt =f) fff p-J?-dv- dt = fy Sff o- E? - dv- dt (18) 


dv 


the distribution of temperature over the entire volume of the conductor is [24]: 


ð z 
YCp =2_y. (KYT) = dy (19) 


where qv is the volumetric energy density of the heat source, ¥ is the specific mass of the conductor, k is the 
thermal conductivity of the conductor, Cp is the specific heat capacity. 

Calculating the temperature distribution on the conductor surface is a very complex procedure because 
the Joule losses depends on the specific resistance of the material, as a result, the surface effect over a constant 
k, also will be depends on the conductor specific resistance. So, the current distribution over the conductor 


Bulletin of Electr Eng & Inf, Vol. 11, No. 5, October 2022: 2510-2517 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 O 2513 


cross section will be depends on the temperature. From these results, we can conclude that the whole process 
of conductor heating is nonlinear and consists of numerous mutual influences [25]. There are some ready-made 
computer packages, such as COMSOL or some other procedures, can be used to solve such complex problems, 
but due to the complexity of this problem, it was decided within this work, for a start, to perform a temperature 
measurement only on the surface of the conductor, which is incomparably simpler than measurements at all 
points of the conductor. To measure the temperature on the conductor surface, a cylinder-shaped temperature 
probe with a diameter of 8 mm and a height of 14 mm has been used. 

The resistance of the probe was measured in connection with four terminals; two voltage and two 
current. The resistance of the temperature probe depending on the temperature is given in Table 1 and 
graphically in Figure 1. FLUKE 9103 dry-well calibrator was used for calibration and maintain a constant 
temperature inside the probe, which changes the temperature in the probe in a range from -25 °C to 140 °C. 
FLUKE 8846a precision multimeter was used to measure the probe resistance. Figure 1 shown that the probe 
has a negative temperature coefficient. 


Table 1. Measured probe resistance 


Nominal temperature 20 40 60 80 100 120 140 160 
Measured probe resistance [kQ] 12.2 5.1 3.5 2.3 1.2 0.75 0.35 0.2 
Measured temperature [°C] 21.3 40.5 60.1 81.3 95 101.9 118.2 132 


Measured probe resistance 


R [K9] 


0 20 40 60 80 100 120 140 
Nominal temperature [°C] 


Figure 1. Dependence of probe resistance on temperature 


2.1. Practical measurement of the conductor surface 

For conducting a measurement of the conductor surface in electrical engineering laboratory, three 
copper conductors with cross-sectional area of: ® = 6 mm?, © = 8 mm? and ® = 12 mm’, were selected, as well 
as one aluminum conductor with cross-sectional area of ® = 8 mm”. Surface temperatures of each conductor 
were measured for 12 different rms values, every 15 seconds, for a total time interval of 30 minutes. When 
measuring the temperature at lower currents, there was no significant heating of the conductor, so the 
measurements could be performed continuously. However, at higher currents, the conductor became more and 
more heated, so before measuring at the next current, the measured conductor had to be cooled by a fan, until 
the conductor temperature reached the same initial value, unique for all measurements. 

By numerical interpolation using the diagrams in Figures 1 and 2, the values of the temperature on 
the surface of the conductor with cross-sectional area (P = 8 mm”) for three different currents are determined 
as shown in Figure 3. The interpolation was done using the MATLAB program, using the "interp1" function. 
The obtained values of individual measurement are shown in Figure 2. Due to the imperfection of the 
measurement procedure, the obtained curves are not smooth as expected. Table 2 illustrates the measuring 
devices used in the practical measurement. 


Table 2. Measuring devices 


No. Measuring devices Purpose 
1 Maxwell MC-25 603A For measuring the effective value of the current. 
2. Power generator for providing simple periodic current up to 180 A. 
3. FLUKE 8846 Multimeter to measure the resistance of the probe. 
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Figure 2. Graph representation of the probe resistance 
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Figure 3. Graph representation of the probe temperature 


3. RESULTS AND DISCUSSION 

Figure 4 illustrates the temperature measurements on the surface of the copper conductor with smallest 
cross section, (P = 6 mm/?), while Figure 5 shows the temperature measurement on the conductor surface of 
the cross-sectional area (® = 8 mm’). In Figures 4 and 5 deviations of the curves at higher currents are observed, 
which is a consequence of the variation of currents due to heating of the conductors and the impossibility of 
compensating for these changes by increasing the supply voltage of the conductors. Namely, the current in the 
conductor is regulated manually, using a potentiometer on the device for generating current, and this manual 
setting could not always monitor changes in temperature and thus changes in current. 
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Figure 4. Graph representation of the surface temperature of the conductor ® = 6 mm? 
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Figure 5. Graph representation of the surface temperature of the conductor ® = 8 mm? 


Figure 6 shows the measurement of the temperature on the surface of copper conductor of the largest cross- 
sectional area (© 12 mm”). In this case, all curves are completely smooth and grow monotonously almost all the 
time, which is a consequence of better maintaining a constant current and a significant increase in the temperature 
of the conductor. The time diagram of the temperature for the aluminum conductor is shown in Figure 7. As shown 
in Figure 7, the decrease in current strength, caused by the increase in specific resistance of the conductor due to the 
increase in temperature, is mostly successfully compensated by a timely increase in the supply voltage of the 
conductor, using a potentiometer on the generator. 
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Figure 6. Graph representation of the surface temperature of the conductor with ® =12 mm? 
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Figure 7. Graph representation of the temperature on the surface of the aluminum conductor 
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4. CONCLUSION 

The aim of the research described in this paper is fully achieved by first determining all the necessary 
characteristics of the unknown temperature probe, and then using it, performed all the planned measurements. 
The measurement results had to be further processed using the procedures of numerical analysis, which was 
also successfully done. The obtained diagrams are very characteristic and can serve as typical forms of heating 
conductors with time-varying currents, including all the effects that accompany such currents. The results 
presented in this paper can also be used to verify the results obtained by calculations. In working on this issue, 
successfully solving all the problems that arose during the work, gained valuable experience in accurate 
temperature measurement. 


REFERENCES 

[I] A. Safdarian, M. Z. Degefa, M. Fotuhi-Firuzabad, and M. Lehtonen, "Benefits of Real-Time Monitoring to Distribution Systems: 
Dynamic Thermal Rating," in JEEE Transactions on Smart Grid, vol. 6, no. 4, pp. 2023-2031, Jul. 2015, doi: 
10.1109/TSG.2015.2393366. 

[2] M. Wang, M. Yang, J. Wang, M. Wang, and X. Han, "Contingency Analysis Considering the Transient Thermal Behavior of 
Overhead Transmission Lines," in JEEE Transactions on Power Systems, vol. 33, no. 5, pp. 4982-4993, Sept. 2018, doi: 
10.1109/TPWRS.2018.2812826. 

[B] X.L. Yu, H. X. Zhang, and M. X. Wang, ‘‘An algorithm for power system state estimation considering line temperature,’’ Zhongguo 
Dianji Gongcheng Xuebao, vol. 38, no. 9, pp. 2561-2570, May 2018, doi: 10.13334/j.0258-8013.pcsee.171161. 

[4] H. Hu, Z. He, X. Li, K. Wang, and S. Gao, "Power-Quality Impact Assessment for High-Speed Railway Associated with High- 
Speed Trains Using Train Timetable—Part I: Methodology and Modeling," in IEEE Transactions on Power Delivery, vol. 31, no. 
2, pp. 693-703, Apr. 2016, doi: 10.1109/TPWRD.2015.2472994. 

[5] L. Battistelli, M. Pagano, and D. Proto, "2$\,\times\,$ 25-kV 50 Hz High-Speed Traction Power System: Short-Circuit Modeling," 
in JEEE Transactions on Power Delivery, vol. 26, no. 3, pp. 1459-1466, Jul. 2011, doi: 10.1109/TPWRD.2010.2100832. 

[6] H. Hu, Y. Shao, L. Tang, J. Ma, Z. He, and S. Gao, "Overview of Harmonic and Resonance in Railway Electrification Systems," in 
IEEE Transactions on Industry Applications, vol. 54, no. 5, pp. 5227-5245, Sept.-Oct. 2018, doi: 10.1109/TIA.2018.2813967. 

[7] J. Q. Zhang and M. L. Wu, ‘‘Power flow algorithm for electric railway traction network based on multiple load port Thévenin 
equivalence,’’ Diangong Jishu Xuebao, vol. 33, no. 11, pp. 2479-2485, Jun. 2018, doi: 10.19595/j.cnki. 1000-6753 .tces. 170501. 

[8] H. Setiadi, A. Swanda, and T. Aryo “Design feedback controller of six pulse three phase rectifier based on differential evolution 
algorithm,” Indonesian Journal of Electrical Engineering and Computer Science, vol. 22, no. 2, pp. 670-677, May 2021, doi: 
10.1159 1Ajeecs.v22.i2.pp670-677. 

[9] S. S. M. Ghoneim, M. Ahmed, and N. A. Sabiha, ‘‘Transient thermal performance of power cable ascertained using finite element 
analysis,” Processes, vol. 9, no. 3, p. 438, 2021, doi: 10.3390/pr9030438. 

[10] C. R. Black and W. A. Chisholm, "Key Considerations for the Selection of Dynamic Thermal Line Rating Systems," in IEEE 

Transactions on Power Delivery, vol. 30, no. 5, pp. 2154-2162, Oct. 2015, doi: 10.1109/TPWRD.2014.2376275. 

[11] F. Wang and X. R. Wang, ‘ʻA power flow analysis method of traction power supply system based on constant—Power Load,”’ 

Proceedings of the CSU-EPSA, vol. 27, no. 3, pp. 59-64, Mar. 2015, doi: 10.3969/j.issn.10038930.2015.03.11. 

[12] H. Lee, C. Lee, G. Jang, and Sae-hyuk Kwon, "Harmonic analysis of the korean high-speed railway using the eight-port 

representation model," in IEEE Transactions on Power Delivery, vol. 21, no. 2, pp. 979-986, Apr. 2006, doi: 

10.1109/TPWRD.2006.870985. 

[13] D. M. Greenwood and P. C. Taylor, "Investigating the Impact of Real-Time Thermal Ratings on Power Network Reliability,” in 

IEEE Transactions on Power Systems, vol. 29, no. 5, pp. 2460-2468, Sept. 2014, doi: 10.1109/TPWRS.2014.2305872. 

[14] J. Ruan, Q. Zhan, L. Tang, and K. Tang, ‘‘Real-Time temperature estimation of Three-Core medium-voltage cable joint based on 

support vector regression,” Energies, vol. 11, no. 6, p. 1405, Mei 2018, doi: 10.3390/en1 1061405. 

[15] K. Chen, Y. Yue, and Y. Tang, ‘‘Research on temperature monitoring method of cable on 10 kV railway power transmission lines 

based on distributed temperature sensor,” Energies, vol. 14, no. 12, p. 3705, Jun. 2021, doi: 10.3390/en14123705. 

[16] L. Xiong, Y. Chen, Y. Jiao, J. Wang, and X. Hu, ‘‘Study on the effect of cable group laying mode on temperature field distribution 

and cable ampacity,’’ Energies, vol. 12, no. 17, p. 3397, Sep. 2019, doi: 10.3390/en12173397. 

[17] J.-A. Jiang, Y. -T. Liang, C. -P. Chen, X. -Y. Zheng, C. -L. Chuang and C. -H. Wang, "On Dispatching Line Ampacities of Power 

Grids Using Weather-Based Conductor Temperature Forecasts," in IEEE Transactions on Smart Grid, vol. 9, no. 1, pp. 406-415, 

Jan. 2018, doi: 10.1109/TSG.2016.2553964. 

[18] I. Theodosoglou, V. Chatziathanasiou, A. Papagiannakis, B. Wie,cek, and G. De Mey, ‘‘El analysis and temperature fluctuations 

*prediction of overhead power lines,’’ International Journal of Electrical Power & Energy Systems, vol. 87, pp. 198-210, May 

2017, doi: 10.1016/j.ijepes.2016.07.002. 

[19] K. J. Ali, G. T. Hasan, and M. A. Ahmed, “Investigate and Analyze the Electromagnetic Field Levels Inside an Electric Power 

Substation,” Tikrit Journal of Engineering Sciences, vol. 24, no. 3, pp. 10-14, Mar. 2017, doi: 10.25130/tjes.24.3.02. 

[20] J. Shu, R. Guan, and L. Wu, ‘‘Optimal power flow in distribution network considering spatial electro-thermal coupling effect,” JET 

Generation, Transmission & Distribution, vol. 11, no. 5, pp. 1162-1169, Mar. 2017, doi: 10.1049/ietgtd.2016.0909. 

[21] K. Kasa8-Lazeti¢é, G. Mijatović, D. Herceg, D. Antić, and M. Prša, ‘‘(in Bosnia: Temperatura na površi provodnika sa vremenski 

promenljivom strujom),” Temperature on the surface of conductors with time-varying current,” 18” International Symposium 

INFOTEH-JAHORINA, pp. 73-77, Mar. 2019. 

[22] L. Kotni, “A proposed algorithm for an overhead transmission line conductor temperature rise calculation,” International 

Transactions on Electrical Energy Systems, vol. 24, no. 4, pp. 578-596, Dec. 2012, doi: 10.1002/etep.1715. 

[23] M. Karahan and O. Kalenderli, “Coupled electrical and termal analysis of power cables using finite element method,” Heat transfer- 

Engeneering Applications, pp. 205-230, 2011. 

[24] X. Dong, C. Kang, Y. Ding and C. Wang, "Estimating the Wind Power Integration Threshold Considering Electro-Thermal 

Coupling of Overhead Transmission Lines," in JEEE Transactions on Power Systems, vol. 34, no. 5, pp. 3349-3358, Sept. 2019, 

doi: 10.1109/TPWRS.2019.2906291. 

[25] S. Frank, J. Sexauer, and S. Mohagheghi, "Temperature-Dependent Power Flow," in IEEE Transactions on Power Systems, vol. 28, 
no. 4, pp. 4007-4018, Nov. 2013, doi: 10.1109/TPWRS.2013.2266409. 


Bulletin of Electr Eng & Inf, Vol. 11, No. 5, October 2022: 2510-2517 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 O 2517 
BIOGRAPHIES OF AUTHORS 


Ali Hlal Mutlaq Oki P he is a lecturer at the Department of electrical and electronic 
engineering, Al-Kitab University, Iraq, where he has been a faculty member since 2005. From 
2015-2018. He graduated with a first-class honours B.Eng. degree in Electronic Engineering 
from Belgrade Universityy, Serbia, in 1982, and M.Sc. in Electrical Engineering 
fromBelgrade University, Serbia in 1984. His research interests are primarily in the area of 
wireless communications and networks as well as in the area of electrical engineering. He 
can be contacted at email: ali33hl@yahoo.com. 


Manhood Ali Ahmed © É P he is an associate professor at the electrical engineering 
Department Al-Shergat engineering college, Tikrit University, Iraq, where he has been a 
faculty member since 2005. He graduated with a first-class honours B.Eng. degree in 
electrical and electronic engineering from Belgrade Universityy, Serbia, in 1984, and M.Sc. 
in electrical engineering from Belgrade University, Serbia in 1986. His research interests are 
in the area of electrical and electronic engineering. He can be contacted at email: 
mh.dham @tu.idu.iq. 


Diadeen Ali Hameed © ÉJ P he is a lecturer in the electrical engineering Department, 
Tikrit university, Tikrit, Iraq. He received the B.Sc. from a assistant professor of electrical 
engineering at University of Tikrit, Iraq. He is an associate professor at the Department of 
electrical engineering, Al-Sherqat engineering college, Tikrit University, Iraq, where he has 
been a faculty member since 2006. He graduated with a first-class honours B.Eng. degree at 
University of Mosul/Iraq and the M.Sc. degrees from University of Yurmouk, Jordin, all in 
computer science and information technology. His research interests are in the area of 
computer and electrical engineering. He can be contacted at email: diaa@tu.edu.iq. 


Ghanim Thiab Hasan Ô ÉJ EJ P he is an associate professor at the Department of electrical 
engineering, Al-Shergat engineering college, Tikrit University, Iraq, where he has been a 
faculty member since 2006. He graduated with a first-class honours B.Eng. degree in 
electrical and electronic engineering from Belgrade Universityy, Serbia, in 1984, and M.Sc. 
in electrical engineering from Belgrade University, Serbia in 1986. His research interests are 
in the area of electrical and electronic engineering. He can be contacted at email: 
ganimdiab @yahoo.com. 


Measurement and analysis of conductor surface temperature independence of current ... (Ali Hlal Mutlaq) 


